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Abstract 

The photoluminescence spectra of band-edge transitions in GaN is studied as a function 
of temperature. The parameters that describe the temperature dependence red-shift of 
the band-edge transition energy and the broadening of emission line are evaluated using 
different models. We find that the semi-empirical relation based on phonon- dispersion related 
spectral function leads to excellent fit to the experimental data. The exciton-phonon coupling 
constants are determined from the analysis of linewidth broadening. 
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In recent years, GaN and its ternary alloys have been extensively studied for their appli- 
cation in blue and ultra-violet light emitting devices, short wavelength lasers, and electronic 
devices for high power and high temperature applications [1]. This tremendous application 
potential originates from the wide range of direct energy band gaps, the high breakdown volt- 
age, and the high saturation electron drift velocity. In spite of technological breakthrough 
in GaN growth, doping and metal-GaN contacting technologies, several fundamental issues, 
such as temperature and doping induced band gap shrinkage(BGS) are still controversial. In 
addition to fundamental interest, these issues are extremely important for the development 
and modeling of GaN-based high temperature and high power devices. The temperature in- 
duced BGS is observed in experiments by a monotonic red shift of either band-to-band(BB) 
and/or excitonic transitions that are observed in bulk as well as low dimensional heterostruc- 
tures. The temperature dependence band gap E g (T) varies from relatively weak in the low 
temperature region to relatively strong at higher temperature region. There is considerable 
controversy regarding the temperature induced BGS and linewidth broadening. All the pre- 
vious studies[2, 3, 4, 5, 6] on GaN have used empirical relations[7, 8] neglecting the role 
of phonon dispersion in GaN[9, 10]to fit the experimental data. There is large variation of 
the values of different BGS and linewidth broadening parameters, for example, in case of 
BGS, the most fundamental parameter a, which is —dE g (T)/dT at high temperature limit, 
varies from 0.4meV/K to 1.2meV/K in the literature[2, 3, 4, 5, 6] and in case of linewidth 
broadening, the most fundamental parameter j L a, which is exciton-longitudinal acoustic 
phonon coupling constant, varies from 0.16meV/K to 0.028meV/K in the literature[4, 6]. 
Passler[ll, 12] has shown the inadequacy of these empirical relations for temperature in- 
duced BGS in semiconductors, in particular the inapplicability[13] of these relations in case 
of wide band gap semiconductors. In addition, residual thermal strain between GaN epi- 
layers and substrates(sapphire and SiC) due to mismatch in lattice constant and thermal 
conductivity is also responsible for large scattering in the values of different BGS parameters. 
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In this letter, we present the temperature dependence of BGS in GaN epilayers from 10K 
to 320K. In particular, the role of background electron concentration on the different param- 
eters responsible for BGS and linewidth broadening is studied. The results are compared 
with different empirical and semi-empirical models for BGS. We have varied the III/V ratio 
to have excess Ga in the buffer layers prior to the growth of top GaN epilayers. Previous 
studies[14, 15] have shown that excess Ga in the buffer layer improves the crystal quality 
of GaN epilayers through stain relaxation and increased adatom surface mobility during 
the initial stage of epitaxial growth. We have recently shown[16] that conductivity and the 
defect-related optical properties of GaN epilayers can be significantly controlled by the III/V 
ratio in GaN buffer layer. 

In this investigation we have used GaN epilayers grown on sapphire substrates by metal 
organic chemical vapor deposition(MOCVD). Prior to the growth of the epilayers, 30nm thick 
GaN buffer layer were grown at 565°C. After growing buffer layer the substrate temperature 
was raised to 1000°C for subsequent 1.6/xm thick GaN epilayers. The background electron 
concentration in the GaN epilayers were varied by different trimethylgallium(TMG) flow 
rate in the buffer layer, while keeping the flow rate of NH 3 and H 2 constant. The TMG, 
NH 3 and H 2 flow rates for all the GaN epilayers were 14/zmol/min, and 1.5 and 2.5slm, 
respectively. The electron concentrations were determined by Hall measurements. We have 
chosen set of undoped n-type samples with background electron concentration in the range 
of 9.8xl0 16 to 2.7xl0 18 cm" 3 . The photoluminescence(PL) spectra were collected in the 
wavelength region of 340-900nm. The samples were kept in a closed cycle He refrigerator 
and were excited with 325nm laser line of He-Cd laser. The PL signal was collected into 
a monochromator and detected with a UV-enhanced Si detector. Fig.l shows the band- 
edge transition, due to the recombination of the exciton bound to neutral donors(known as 
I 2 exciton) [17, 18], which can be observed at higher temperature(>300K) due to relatively 
large binding energy(~25meV). To avoid the ambiguity in determination of band gap due 
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to transition from excitonic to other band-edge transitions at higher temperature, we have 
limited our experiment till 320K. 

Temperature dependence of the PL spectra of two n-type GaN samples with different 
background electron concentrations is shown in Fig.l, which shows a maximum redshift of 
about 65meV as temperature increases from 10K to 320K. An appropriate fitting function is 
required to obtain material-specific parameters from experimentally measured E g (T). There 
are two set of fitting function for E g (T) in the literature: (i) empirical relations proposed by 
Varshni[7] and Vina et al[8] and (ii) semi-empirical relations based on the electron-phonon 
spectral function /(e) and the phonon occupation number n(T), e is the phonon energy . 

The most frequently used empirical relation for numerical fittings of E 9 (T) was first 
suggested by Varshini[7] and given by 

E g {T) = EM - -0^ (1) 

where E 9 (0) is the band gap at OK, a is the T — > oo limiting value of the BGS coefficient 
dE g (T) I dT and (3 is a material specific parameter. This model represents a combination 
of a linear high-temperature dependence with a quadratic low-temperature asymptote for 
E g (T). Though this phenomenological model gave reasonable fittings of E 9 (T) in elemental, 
III-V and II- VI semiconductors with E g < 2.5eV, there have been several problems with 
this relation, for example, (i) this relation gives negative values of a and (5 in case of wide 
band gap semiconductors [7, 13], (ii) j3 is a physically undefinable parameter believed to 
be related to Debye temperature Od of the semiconductor, but this connection has been 
doubted strongly in several cases [12] and (iii) it has been shown that this relation cannot 
describe the experimental data of E 9 (T) even in GaAs[19]. Fig. 2 (a) shows the comparison of 
experimental E 9 (T) with Varshni's relation[7]. We have obtained the values of a in the range 
of 0.54 meV/K to 0.63 meV/K and (3 in the range of 700 K to 745K for the samples with 
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electron concentrations of 9.8xl0 16 to 2.7xl0 18 cm" 3 . The parameter (3 in this model only 
gives an estimation of Od, which is about 870K in GaN. It is clear that though the fitting is 
good in the low temperature region(< 100K), but fitting is poor in the intermediate(~100K) 
and high temperature region(> 200K). 

Vina et al[8] first emphasized that total BGS AE g (T) = E g (0) — E g (T) is proportional 
to average phonon occupation numbers n(T) = exp{-^f) — 1 and proposed an empirical 
relation, which can be expressed as 

E ' iT) = EM ~ ^fpi (2) 

where cub = ff^-, is a material parameters related to electron- phonon interaction, Q B = 
hul k e * f represents some effective phonon temperature and uj e ff is the effective phonon fre- 
quency. Though this phenomenological model gave reasonable fittings of E S (T) in different 
semiconductors, there are also several problems with this relation, for example, (i) at low 
temperature, this model shows a plateau behavior of E g (T), which is not observed exper- 
imentally, (ii) at higher temperature(>50K), this model predicts AE g (T) oc exp(—Q B /T), 
but, in most cases AE g (T) oc T 2 is observed experimentally, (iii) it has been shown that this 
relation cannot describe the experimental data of E S (T) even in GaAs[19]. Fig. 2(b) shows 
the comparison of experimental E 9 (T) with Vina's relation. We have obtained the value of 
B in the range of 400K to 450K and a B in the range of 0.40meV/K to 0.44meV/K for 
the samples with electron concentrations of 9.8xl0 16 to 2.7xl0 18 cm" 3 . There is no doubt 
that this gives better fitting than Varshini's relation in the high temperature region, but the 
fitting in the low temperature region is poor. 

As mentioned earlier, it can be shown that the contribution of individual phonon mode 
to the temperature induced BGS is related to average phonon occupation number n(T) and 
electron-phonon spectral function /(e). Essentially, E g (T) can be analytically derived from 
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the expression 



E g (T) = E g (0)- Jdef(e)n(e,T) (3) 

The electron-phonon spectral function /(e) is not known a priori and extremely complicated 
to calculate from the first principle. The other option is to use different approximate function 
for /(e) to derive temperature dependence of BGS. It has been emphasized conclusively [11, 
12] that the indispensable prerequisite of estimation of different parameters obtained from 
the experimentally measured E g (T) is the application of an analytical model that accounts 
for phonon energy dispersion. The band gap shrinkage results from the superposition of 
contributions made by phonons with largely different energies, beginning from the zero- 
energy limit for acoustical phonons up to the cut-off energy for the optical phonons. The 
basic features of phonon dispersion 5 p h and the relative weight of their contributions to E g (T), 
may vary significantly from one material to the other. The curvature of the nonlinear part 
of the E g (T) is closely related to the actual position of the center of gravity, e(e = hu), 
and the effective width Ae, of the relevant spectrum of phonon modes that make substantial 
contribution to E g (T) and this has been quantified by phonon-dispersion co-efficient 5 p h{= 
4^). It has been shown[13] that above two empirical models(Eqn.l and Eqn.3) represent the 
limiting regimes of either extremely large(5 p / l ~1) in case of Varshini's relation or extremely 
small(<5 p?i ^0) in case of Vina's relation. Both these models contradict physical reality 
for most semiconductors, whose phonon dispersion co-efficient vary between 0.3 to 0.6[13]. 
Several analytical models have been presented using different forms of the spectral function 
/(e). Passler[ll] has proposed the most successful model, which takes power-law type spectral 
function, /(e) = Vj^ and the cut-off energy eo, which is given by eo = ^^-ksQp- 

Here, v represents an empirical exponent whose magnitude can be estimated by fitting the 
experimental data E g (T). Inserting, the spectral function /(e) and cut-off energy e into the 
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general equation for band gap shrinkage (Eqn. 2), Passler[ll] obtained an analytical expression 
for E g (T), given by 



where p=z/+l and a p is the oo limit of the slope dE ff (T)/dT, Q p is comparable with the 
average phonon temperature [11], 6 P ~ e/k B , and the exponent p is related to the material- 
specific phonon dispersion co-efficient S p h, by the relation 5 p h ~ 1/Vp 2 ~~ 1- Depending on 
the value of 5 p h, there are regimes of large and small dispersion which are approximately 
represented within this model by exponents p < 2 and p > 3.3 respectively. Fig.2(c) shows 
the comparison of experimental E g (T) with Eqn. 4. We obtain an excellent fit to the experi- 
mental data. The value of a p obtained here is in the range of 0.50meV/K to 0.54meV/K and 
the value of O p in the range of 500 K to 510 K for the samples with electron concentrations 
of 9.8xl0 16 to 2.7xl0 18 cm -3 . The value of the exponent p which is related to the material 
specific degree of phonon dispersion is obtained between 2.50 to 2.65 for different GaN sam- 
ples. The value of p obtained here lies in the intermediate dispersion regime and the value 
of the ratio Q p /Qd is 0.58, which means that the center of gravity of the relevant electron- 
phonon interaction are located within the upperhalf of the relevant phonon spectra[9, 10]. 
Fig. 3 shows that the total redshift AE g = E g (10K) — E g (320K) and the phonon dispersion 
co-efficient 5 p h decrease with electron concentration. This trend may be explained in terms 
of the screening of the electron-phonon interaction with higher electron concentration in the 
heavily doped samples. 

Fig. 4 shows the temperature dependence of full width at half maximum(FWHM) of 
emission PL line in GaN. The FWHM, T(T) can be given by[20] 






(5) 
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where r is the intrinsic linewidth at T=OK due to temperature-independent mechanisms, 
such as defects and electron-electron interaction, and T LO is the strength of the exciton-LO- 
phonon coupling and Ti^lo is the LO phonon energy. Fig.4 shows the fit to the experimental 
data using Eqn.5. Insets of Fig.4 shows how 7^ and T LO vary with electron concentra- 
tions. The estimated LO-phonon energy(~92meV) is very close to the value determined by 
independent measurements[10]. We find exciton-phonon coupling constants are larger than 
other III-V semiconductors. It has been found that both exciton-LA phonon and exciton- 
LO phonon interactions become weaker with increasing electron concentration, as shown in 
Fig.4, signifying the important role of exciton-phonon or electron-phonon interaction on the 
overall redshift due to temperature induced BGS. These findings corroborate with results 
shown in Fig. 3. 

In conclusion, we have measured temperature induced band gap reduction and broadening 
of linewidth in GaN using photoluminescence spectroscopy. The importance of electron- 
phonon interaction on the band gap shrinkage has been established. It has been found that 
phonon-dispersion based semi-empirical relation is required to explain the experimental data. 
Screening of electron-phonon interaction is responsible for decrease of temperature induced 
BGS in samples with higher electron concentration. 
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Figure Captions 



Figure 1. PL spectra of band-edge excitonic transition at different temperatures from 
fOK to 320K in GaN epilayers with electron concentration of (a) 1.3xl0 18 cm~ 3 and (b) 
3.8xl0 17 cm- 3 . 

Figure 2. Temperature dependence of the peak positions of band-edge excitonic transition 
in GaN sample with electron concentration of 3.8x 10 17 cm~ 3 . Similar fits have been ob- 
tained in case other samples with electron concentration of 9.8x 10 16 cm -3 , 1.3xl0 18 cm~ 3 
and 2.7x 10 18 cm~ 3 . Solid lines are fit to experimental data with (a) Eqn.l, (b) Eqn.3 and 
(c) Eqn.4. 

Figure 3. Electron concentration dependence of BGS, AE g = E g (10K) — E g (320K). Inset 
shows how the phonon-dispersion parameter S p h derived from relation 5 p h = 1 / V 1 — p 2 
changes with electron concentration n. Solid lines are guide for eyes. 

Figure 4. Temperature dependence of the FWHM of PL emission line for GaN epilayer 
with electron concentration of 3.8xl0 17 cm~ 3 . Solid line is the fit to experimental data with 
Eqn.5. Similar fits have been obtained in case other samples with electron concentration of 
9.8xl0 16 cm~ 3 , 1.3xl0 18 cm~ 3 and 2.7xl0 18 cm~ 3 . Upper inset shows how electron-acoustic 
phonon coupling constant ('J la) varies with electron concentration n. Lower inset shows how 
electron-optical phonon coupling constant(rLo) varies with electron concentration n. Solid 
lines are guide for eyes. 
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